ApJ in press 

Preprint typeset using LAT^X style cmulateapj v. 5/2/11 



THE OPTICAL AND NEAR-INFRARED TRANSMISSION SPECTRUM OF THE SUPER-EARTH GJ1214B: 
FURTHER EVIDENCE FOR A METAL-RICH ATMOSPHERE 

Jacob L. Bean 1,8 , Jean-Michel Desert 1 , Petr Kabath 2 , Brian Stalder 1 , Sara Seager 3 , Eliza Miller-Ricci 
Kempton 4,8 , Zachory K. Berta 1 , Derek Homeier 5 , Shane Walsh 6 , & Andreas Seifahrt' 

ApJ in press 



o 

» 

o 
O 

(N 

w 

6 

CO 

03 



> 

oo 
in 
o 

ON 

o 



13 



ABSTRACT 

We present an investigation of the transmission spectrum of the 6.5 M® planet GJ 1214b based 
on new ground-based observations of transits of the planet in the optical and near-infrared, and 
on previously published data. Observations with the VLT + FORS and Magellan + MMIRS using 
the technique of multi-object spectroscopy with wide slits yielded new measurements of the planet's 
transmission spectrum from 0.61 to 0.85 /an, and in the J, H, and K atmospheric windows. We 
also present a new measurement based on narrow-band photometry centered at 2.09 /um with the 
VLT + HAWKI. We combined these data with results from a re-analysis of previously published FORS 
data from 0.78 to 1.00 urn using an improved data reduction algorithm, and previously reported values 
based on Spitzer data at 3.6 and 4.5 /im. All of the data are consistent with a featureless transmission 
spectrum for the planet. Our iT-band data are inconsistent with the detection of spectral features at 
these wavelengths reported by Croll and collaborators at the level of 4.1 a. The planet's atmosphere 
must either have at least 70% H2O by mass or optically thick high-altitude clouds or haze to be 
consistent with the data. 

Subject headings: planets and satellites: atmospheres — planets and satellites: individual: GJ 1214b 
— techniques: photometric 



1. INTRODUCTION 

Recent results from the Kepler mission indicate that 
exoplanets intermediate in size between the terrestrial 
and the ice giant planets in the Solar System are common 
at small orbital separations from main sequence stars 
(jBorucki et all 120111: IBoward et all 120111) . However, de- 
spite the wealth of knowledge of the sizes and masses 
for these planets that will be obtained from the Kepler 
data and follow-up efforts, there will still be a substan- 
tial gap in our understanding of these objects. This is 
because knowledge of only the mass and radius of plan- 
ets in this regime is not sufficient to determine their bulk 
compositio n due to degenerac i es that exist in theoret i- 
cal models ([Adams et al.ll20~0 8: Rogers fc Seagenl27LT0al ). 
Specifically, most planets with masses and radii in the 
intermediate-size range can be modeled equally well by 
substantially different combinations of mass and compo- 
sition for the interior and surrounding atmosphere. 

Fortunately, there is a way to break the degen- 
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eracics between different models for a planet when 
a unique solution for the bulk composition can not 
be found. The method is to obta in constraints on 
the p l anet's atmospheric comp osition ([Miller-Ricci et all 
120091 : iMirler-Ricci fe Fortnevi 120101 ). This information 
provides a boundary condition to further constrain both 
rigorous physical models and also cosmogonical argu- 
ments about likely origin and evolutionary scenarios. 
Therefore, detailed characterization of representative ob- 
jects to further constrain their bulk compositions is an 
important complement to the statistical study of the oc- 
currence rate of populations that is being enabled by 
Kepler. 

Th e transiting planet GJ 1214b (Cha rbonneau et al.l 
120091 ) is one object that is potentially a key to unlocking 
the mystery of the bulk compositions of intermediate- 
size planets using atmospheric studies. Its determined 
mass (6.5 M®) and radius (2.65 R®) put it firmly in the 
degenerate region of parameter space for interior and at- 
mosphere mass and compositions. The only thing that 
can be said with a high degree of confidence based on 
just the mass and radius of the planet is that it must 
have a substantial atmosphere because it is too large to 
be composed solely of solid m aterial ([Rogers fe Seagerl 
2010 d iNettelmann et alll2011l) . 

Ro gers fc Seagerl (|2010bD identified three possible ori- 
gins for the gas layer on GJ 1214b and they presented dis- 
tinct interior structure models based on these scenarios 
that are consistent with its measured mass and radius. 
The model scenarios are (1) a "Mini-Neptune" planet 
that is composed mainly of solid rock and ice and that has 
a significant primordial atmosphere accreted from the 
protoplanetary nebula and with a composition roughly 
similar to that of the Sun (i.e., hydrogen-dominated), 
(2) a "Water World" planet composed primarily of water 
ice and having a secondary gas envelope formed by sub- 
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limation and dominated by water vapor, and (3) a true 
"Super-Earth" planet composed of purely rocky material 
with a secondary atmosphere formed by outgassing and 
probably composed mainly of hydrogen. Distinguish- 
ing between these competing models for GJ 1214b would 
have important consequences because they each imply 
a very different formation and evolutionary history for 
the planet. Constraining the formation and evolution- 
ary history for this planet would be an important step 
in our quest to obtain a general understanding of the 
intermediate-size planets for which this object appears 
to be an archetype. 

The most interesting characteristic of GJ 1214b is that 
it orbits a very small M dwarf, and thus the system has 
a planet-to-star radius ratio comparable to a Jupiter-size 
planet transiting a Sun-like star. This quality makes it 
the most feasible known intermediate-size planet for at- 
mospheric studies using transit spectroscopy techniques. 
The three models proposed for the planet would also 
exhibit very different transmission spectrum features. 
Mini-Neptune and true Super-Earth planets with their 
hydrogen-dominated, and thus large scale height, atmo- 
spheres would exhibit relatively large spectral features in 
transmission owing to absorption by trace gases like wa- 
ter and methane, and scattering by molecular hydrogen 
and cloud or haze particles. On the other hand, a Water 
World planet with a primarily water vapor atmosphere, 
and thus small scale height, would exhibit relatively small 
spectral features in transmission. 

We recently obtained measurements of GJ 1214b's 
transmission spectrum in the red optical (0.78 - 1.00 /im) 
with the FORS instrument on the VLT using a new 
ground-based techniq ue to constrain th e planet's atmo- 
spheric composition (Bean et all 120101 ) . This remark- 
ably featureless spectrum is consistent with a model 
for a metal-rich atmosphere with a small scale height, 
and inconsistent with a model for cloud-free hydrogen- 
dominated atmospheres, which would have a large scale 
height. Clouds or haze in a hydrogen-dominated at- 
mosphere could conceivably also yield a flat spectrum 
consistent with the observations. Subsequently, we pre- 
sented broadband photometric measurements obtained 
at 3.6 and 4.5 /im with Spitzer that provided addi- 
tional constraints o n the composition of the atmosphere 
(|Desert et al.l 120 111 ) . These data were also consistent 
with a featureless transmission spectrum when analyzed 
in isolation and in combination with the FORS data, but 
the overall interpretation of the flat spectrum remained 
uncertain due to the possibility of clouds or haze pro- 
viding a gray opacity source, and non-equilibrium abun- 
dances of methane, which is the primary expected opac- 
ity source at 3.6 /im. 

In cont rast to the result s showing a featureless 
spectrum, ICroll et all ()2011aD presented ground-based 
near-infrared photometry observations that indicated 
GJ1214b's transmission spectrum has a large feature 
in the _R" s -band (bandpass approximately 2.0 - 2.3 /im) 
based on the measurement of a substantially deeper tran- 
sit in that band as compared to the observations at other 
wavelengths. The main interpretation of the observations 
in this case is unambiguous; the large scale height of a 
hydrogen-dominated atmosphere is required for such a 
fe ature to be obse r ved. 

iCrossfield et a l. (2011) presented high-resolution spec- 



troscopy of GJ 1214b between 2.1 - 2.4 /im. They de- 
tected no spectral features, but could not comment 
on the absolute depth of the transit at these wave- 
lengths due to the pure ly differential nature of their data. 
I Crossfield eTall (l20TlV came to similar conclusions as 
IDesert et al.l f|201 lh did based on the Spitzer data be- 
cause methane is also the expected source of wavelength- 
dependent opacity in the their observational window. 

The observational results showing a featureless trans- 
mission sp e ctrum for GJ 1214b and the result from 
ICroll et "all (|2011aD indicating large spectral features in 
the near-infrared are not strictly in conflict because 
the y either wer e obta ined at different w a velen gths [i.e., 
the iBean et al.l (|2010l ) and IDesert et al] (120111) studies! 
or ha ve different sensitivities [i.e.. the ICrossfield et al.l 
(| 2 1 If) investigation]. Also, there is a plausible qualita- 
tive explanation for all the observations: a hyrodgen- 
dominated atmosphere depleted in methane and with 
clouds or haze opaque at optical wavelengths. However, 
there is some tension because, for all the plausibility of 
the qualitative model, it is not supported by physical cal- 
culations at this point, and it has to be finely tuned to 
match the observations. 

Given the outstanding issues in our understanding of 
GJ 1214b's atmosphere, we were motivated to conduct 
further observational studies of its transmission spectrum 
to improve our understanding of this important world. In 
particular, we aimed to make independent measurements 
of the planet's transmission spectru m in the if -band re- 
gion using a different technique than ICroll et al.l (|2011aD 
used, and also to make measurements further in the blue 
part of the optical to search for the signature of scatter- 
ing from cloud, haze, or gas particles. We present here 
the results of these investigations. In §2 we describe our 
observations and data reduction. In §3 we describe our 
analysis of the data. We compare our new measurements 
of the planet's transmission spectrum to previous mea- 
surements and theoretical models in §4. We conclude in 
§5 with a summary of the results. 

2. DATA 
2.1. Magellan + MMIRS 
2.1.1. Observations 

We observed transits of GJ 1214b on 2011 May 15 and 
18 using the MMIRS instrument (jMcLeod et al.l 120041) 
on the Magellan (Clay) telescope at Las Campanas Ob- 
servatory. We used the same multi-object spectroscopy 
approach for the MMIRS observations as was used for 
the FORS observations of G J 1214b that we presented 
previously ()Bean et al.ll2010D . We gathered time-series 
spectra for GJ 1214b and three other reference stars of 
similar brightness within 6' using slits with lengths of 30" 
and widths of 12". The off-axis guide and wavefront sen- 
sor of the instrument, which is sensitive to wavelengths 
between 0.6 and 0.9 /im, was used to feed corrections to 
the telescope control and active optics systems. We used 
an HK grism as the dispersive element and an HK filter 
to isolate the first order spectra. Complete spectra from 
1.23 to 2.48 /im with a dispersion of 6.6Apixel _1 were 
obtained for all the objects. 

As described below, only the data obtained on the sec- 
ond night are reliable. The observations for this transit 
spanned 2.63 h from UT 06:15 to 08:53. Exposure times 
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Fig. 1. — An example spectrum extracted from the MMIRS data 
for GJ 1214. The broad bandpasses in the J, H, and K atmospheric 
windowns are indicated by the grey boxes. 

were 25 s, and the overhead per exposure including the 
read and reset time of the detector was 14 s. A total of 
239 exposures were obtained, 79 of which were in tran- 
sit. The observations began when the field was at an 
airmass of 1.20, and by the end the field had set to an 
airmass of 1.58. Conditions were clear and the seeing 
varied between 0.5 and 0.9" as estimated by the width 
of the spatial profiles of the obtained spectra. 

2.1.2. Data reduction 

The MMIRS data were recorded using the "up-the- 
ramp" sampling mode of the detector in which non- 
destructive reads were obtained every 5 s and stored in 
a data cube. In the first step of the data reduction, we 
subtracted the first of these reads from the others on 
a pixel-by-pixel basis in each frame to remove the bias 
and persistence from previous exposures. The estimated 
dark current of the detector is 0.06 e - s -1 pixel -1 , which 
is negligible for this experiment. After subtracting the 
first read in the data cube for an exposure, we then fit a 
slope to the up-the-ramp samples for each pixel to deter- 
mine the count rate. We arrived at the final estimate of 
the total pixel counts by multiplying the count rate by 
the total exposure time. 

The reduction of the MMIRS data cubes gave reason- 
able results for the second data set, but we noticed a re- 
peating pattern of discontinuities in the residuals when 
fitting the up-the-ramp samples for the first data set. 
The origin of this effect is likely a problem in the detec- 
tor electronics. We chose not to include the first data 
set in the final analysis because of this problem; in what 
follows we focus only on the analysis of the data set ob- 
tained on the 18th. 

After the data processing to collapse the data cubes, 
we applied spectroscopic flat-field corrections, subtracted 
the background, and extracted one-dimensional spectra 
for the objects using the same algorithms as f or our previ- 
ous m ulti-object spectroscopy observations (iBean et "all 
2010). We established the wavelength scale of the ex- 
tracted stellar spectra based on spectra of an Ar emission 
lamp obtained following the transit observation using a 
copy of the science mask with the slit widths set to 0.5". 
We adopted a third-order polynomial form for the wave- 
length solution. The wavelength calibration was done 
independently for each of the four slits, and the maxi- 
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Fig. 2.— The MMIRS Tf-band raw time-series data for GJ1214 
(filled circles) and the composite reference (open diamonds). 

mum dispersion from the fits was 0.3 A. After applying 
the wavelength solution to the data, we noticed offsets 
between the positions of telluric absorption lines in the 
spectra extracted from the first image. These offsets are 
consistent with a misalignment in the rotation angle be- 
tween the science mask and the calibration mask. We de- 
termined and applied offset corrections to the zero point 
of the wavelength scale of each object. The maximal 
correction size was 38 A, which is equivalent to approxi- 
mately six pixels. 

The positions of the spectra on the detector varied in 
both the spatial and spectral directions during the time- 
series due to imperfect guiding and instrument flexure, 
and this led to apparent shifts in the extracted Id spec- 
tra relative to each other. We cross-correlated the spec- 
tra for each object relative to the first spectrum in the 
time-series, and then shifted them by this amount to cor- 
rect for the effect of the drift. The position of the spec- 
tra drifted by 2.5 pixels over the course of the first 25 
exposures, and then there was a sudden shift of approx- 
imately 3 pixels. The position of data after the large 
jump remained relatively stable with a rms of 0.15 pix- 
els and only a small long-term drift totaling 0.13 pixels. 
We removed the first 30 points from the time-series to 
avoid problems caused by the large instability at the be- 
ginning. There is still sufficient out-of-transit baseline 
(45 minutes) with these points removed, and the model 
fits are substantially better for the trimmed light curves. 

An example spectrum of GJ 1214 is shown in Figure[TJ 
The signal-to-noise ratio (SNR) in the GJ 1214 spectra 
estimated from photon counting statistics reaches up- 
wards of 475 pixel -1 in the if-band atmospheric window 
where the maximum number of counts were recorded. 
We estimate that the resolution of the stellar spectra is 
approximately 50 A (7.5 pixels). This estimate is based 
on the intrinsic resolving power of the grism (R = A/AA 
?» 560 for a 1.0" slit), the atmospheric seeing, errors in 
the correction for the motion of the stars on the slits dur- 
ing the observations, and the uncertainty in the absolute 
wavelength scale. 

2.1.3. Extraction of the photometric time-series 
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Fig. 3. — Normalized residuals (circles) from the MMIRS broadband i-sT-band light curve fits without the airmass (left panel) and the 
spatial spectrum position (right panel) decorrelations as a function of those parameters. The decorrclation fits are indicated by the solid 
lines. 



We created broadband photometric time-series from 
the extracted spectra in the J, H, and K atmospheric 
windows, and four spectroscopic time-series in the K- 
band of width 0.1 /jm by summing the obtained spectra 
within a bandpass. As discussed in §3, both the broad- 
band and spectroscopic data provide insight to the com- 
position of the planet's atmosphere. We do not feel that 
the J- and ii-band data are of sufficient quality to justify 
subdividing them into spectroscopic channels. 

The definitions of the utilized bands and the deter- 
mined planet-to-star radius ratios (i? p /i?*) from the light 
curve modeling (see §3) are given in Table [TJ The limits 
of the broad bandpasses were chosen to match standard 
near-infrared filter curvetQ as closely as possible to avoid 
regions with significant telluric water absorption and to 
enable quick comparison with oth er results, i n partic - 
ular the if s -band measurement of iCroll et al.l (|2011aD . 
However, our data do not include the entire range of the 
standard J-band, and the combination of the grism and 
order-blocking filter we used has substantially different 
transmission properties than broadband filters. The ex- 
act locations of the bandpass edges were also set so that 
the spectra were not cut in the middle of strong absorp- 
tion lines. Therefore, appropriately weighted integration 
of the atmospheric models for the planet are necessary for 
robust comparison of results even at similar wavelengths. 

The time-series data for GJ 1214 in each bandpass were 
corrected by dividing out the sum of the data for the ref- 
erence stars. In the subsequent light curve modeling, 
we found that including all of the reference stars in this 
step yielded a worse correction than when one of them 
was excluded. A possible explanation for this is that 
the spectra for the problematic reference star fell along 
the boundary between different read-out channels in the 
detector. The motion of the spectra during the time- 
series resulted in the spectra for this object moving back 
and forth across the boundary and this could have intro- 
duced additional noise in the data. For the final analysis, 



1 See http: //www. cf ht .hawaii . edu/Instruments/Filters/wircam 

for the properties of the filters used by (Jroll ct al. (20fla). 



TABLE 1 

Photometric Bandpasses & Transit Depths for the 
MMIRS Data 



Band 


Wavelength (;im) 


Rp / R-k a 


Broadband analysis 






J 


1.25 - 1.33 


0.1158 ±0.0024 


H 


1.50 - 1.77 


0.1146 ±0.0014 


K 


1.98 - 2.31 


0.1158 ±0.0006 


Spectroscopic analysis 






channel 1 


1.98 - 2.08 


0.1156 ±0.0007 


channel 2 


2.08 - 2.18 


0.1163 ±0.0007 


channel 3 


2.18 - 2.28 


0.1158 ±0.0006 


channel 4 


2.28 - 2.38 


0.1163 ±0.0011 


a From an analysis with the transit parameters 


a I R* and i fixed 



to 14.97 and 88.94°, respectively. 

we excluded the data for the bad reference star, and we 
used only the data for the other two reference stars to 
correct the GJ 1214 time-series. The estimated photon- 
limited uncertainties in the photometry for the reference 
star were propagated through to the final uncertainties 
in the GJ 1214 data. An illustration of the reference star 
correction to the broadband if-band data is shown in 
Figured 

After correcting the time-series for GJ 1214 using the 
reference star data, we identified additional systematics 
correlated with the airmass of the observations and the 
positions of the spectra in the spatial dimension on the 
detector. We fit for linear decorrelations against these 
two parameters simultaneously with the light curve mod- 
eling described in §3 to remove these effects (two free 
parameters in addition to the normalization). Adding 
higher-order terms and decorrelating against airmass and 
spatial pixel position for the individual stellar fluxes (i.e., 
GJ 1214 and the reference stars) did not improve the 
results. The correlations between the reference star- 
corrected light curve residuals and the resulting linear 
fits for these two parameters are shown in Figure [3] 
htiT| The MMIRS broadband light curves for GJ 1214 after 
all corrections and decorrelations exhibit residuals from 
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Fig. 4. — An example spectrum extracted from the FORS blue 
data for GJ 1214 

a best-fit model with rms of 1532, 887, and 706 ppm in 
the J-, H-, and JT-bands, respectively. These values are 
factors of 3.9, 6.0, and 3.4 higher than expected from the 
photon-limited uncertainties. The if-band data exhibit 
an especially larger-than-expected scatter of unknown 
origin. We experimented with changing the bandpasses 
to investigate this issue. Shifting the bandpasses and also 
using smaller and larger windows did not significantly 
affect the results. We speculate that perhaps the trans- 
parency of the Earth's atmosphere varies on different spa- 
tial scales among the near-infrared atmospheric windows; 
and that this scale, along with the number, color, and 
spatial distribution of reference stars, sets the fundamen- 
tal limit when working in the < 1000 ppm photon-limited 
regime in the near-infrared. When binned over 60s, 
the rms of the if-band residuals is 443 ppm. These are 
the highest-precision ground-based near-infrared transit 
light curve data that we are aware of. The rms of the 
residuals for the iC-band spectroscopic channels range 
from 923 ppm to 1552 ppm, and the scaling from the ex- 
pected levels of noise range from 2.3 to 2.8. 

2.2. VLT + FORS 

We observed one transit of GJ 1214b on 2011 July 3 
using the m ulti-object spectroscopy mode of the FORS 
instrument (jAppenzeller et al.lll998l ) on UT1 of the VLT 
facility at Paranal Observatory. The instrument setup 
and data reduction were very similar to the pre vious 
FORS observation s described in IBean et al.l ((20100 and 
IBerta et all (|2011D . and we limit our discussion here to 
the unique aspects of the current study. 

We used the 600RI grism with the GG435 filter to ob- 
tain spectra over the range 610 - 850 nm for GJ 1214 and 
five reference stars. Observations began at UT 02:04 and 
continued for 3.12 h to UT 05:11. The field was at an air- 
mass of 1.20 at the beginning, rose to an airmass of 1.15, 
and then set to an airmass of 1.33 by the end. The ex- 
posure times were 18 s, which gave a total duty cycle of 
55 s including 37 s of overhead. A total of 194 exposures 
were obtained, 57 of which were during the transit. 

Our approach for the data reduction, spectral extrac- 
tion, and wavelength calibratio n were the same as for 
the previous FORS observations (|Bean et al.| [2010). with 
one exception. During the wavelength calibration step, 
we identified a few lamp lines in the FORS atlas that 
must have inaccurate assigned wavelengths because they 
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Fig. 5. — The broadband FORS blue raw time-series data for 
GJ 1214 (filled circles) and the composite reference (open dia- 
monds). 

were significantly discrepant from a low-order polyno- 
mial dispersion model. When we did not include these 
lines in the solution, it became clear that a third-order 
polynomial model for the dispersion was justified instead 
of the linear model we had used previously and that is 
also advocated in the instrument manual. The residuals 
from the third-order dispersion fit with the bad lines ig- 
nored are typically 0.005 nm, whereas the residuals are 
typically 0.5 nm from a linear model with the bad lines 
included. 

This new insight to the wavelength scale of FORS 
grism data prompted us to re-analyz e the previous dat a 
for the two transits of GJ 1214b from IBean et al.l ([2010), 
and we present the new results below. We believe that 
the resolution of the data can be reliably increased now 
that the wavelength solution has been significantly im- 
proved, and we re-bin the old spectroscopic data to 10 nm 
instead of the 20 nm used before. The reduction and 
correction of the data are otherwise unchanged. The 
basic properties of the data remain the same af t er the 
re-analysis, and the conclusions from IBean et al.l (|2010f ) 
are strengthened. In the rest of the paper we refer to 
the new FORS data as the FORS "blue" data, and the 
re-analyzed old data as the FORS "red" data. 

An example spectrum of GJ 1214 extracted from the 
FORS blue data is shown in Figure The signal-to- 
noise ratio in the GJ 1214 spectra estimated from pho- 
ton counting statistics reaches upwards of 300pixel _1 at 
900nm, but only 50pixel _1 at 600nm due to GJ 1214's 
extremely red spectral energy distribution. 

We created both a broadband and spectroscopic light 
curves from the FORS blue data by summing over the full 
wavelength range and over twelve channels with widths 
of 20 nm, which is the lowest we feel that we can go 
without systematic errors given the quality of the data. 
We combined the data for four of the reference stars and 
divided the composite from the GJ 1214 time-series to 
correct it. The fifth reference star was found to yield a 
worse correction to the photometry of GJ 1214 and was 
not included in our final analysis. An illustration of the 
reference star correction to the broadband data is shown 
in Figure [5] 
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Fig. 6.— The HAWKI raw time-series data for GJ 1214 (filled 
circles) and the composite reference (open diamonds). 

As with our previous study of FORS data, we identified 
a slow and smooth variation in the corrected relative flux 
values with time superimposed on the expected transit 
light curve morphology. This effect is probably due to the 
smudges on the FORS atmospher ic dispersion corrector 
described bv lMoehler et al.l (|2010l) . We fit for a second- 
order polynomial as a function of time simultaneously 
with the transit modeling (see §3) to account for this 
(two free parameters in addition to the normalization). 

The FORS blue broadband light curve for GJ 1214 af- 
ter all corrections and decorrelations exhibits residuals 
from a best- fit model with a rms of 457 ppm. This is 
a factor of 2.5 larger than expected from the estimated 
photon noise. The rms of the residuals for the spectro- 
scopic channels range from 665 ppm to 1785 ppm, and 
the scaling from the expected levels of noise range from 
1.2 to 1.4. 

2.3. VLT + HAWKI 

We observed one transit of GJ 1214b on 2011 August 10 
using the standard imaging mo de of the HAWKI instru- 
ment (|Kissler-Patig et alJl200SD on UT4 of the VLT facil- 
ity at Paranal Observatory. The observations were done 
using the NB2090 filter (A c = 2095nm, width = 20 nm). 
In contrast to previously reported HAWKI observa- 
tions of exoplanet t r ansits (e.g., iGillon et al.l 120091 : 
lAnderson etaLl [2010t iGibson et ail [2010). we chose to 
maintain a single pointing throughout the observations 
instead of using a dither pattern. This "staring" ap- 
proach has been used for observations with a similar in- 
strument to obtain the highest-precision ground-based 
transit light curves in the near-infrared up to now (a p- 
proxim ately 700 ppm minute -1 ; e.g.. lCroll et al.ll2~010al lb1. 
2011b). We reach a similar level of precision in our data 
using the staring approach, and this suggests that the 
single pointing method is indeed the best technique to 
use for exoplanet transit observations with HAWKI. 

Observations began at UT 01:03 and continued for 
2.68 h to UT 03:44. The field was at an airmass of 1.15 at 
the beginning, and it set to an airmass of 1.64 by the end. 
The exposure times were 1.6762 s, which is the minimum 
possible for the normal mode of HAWKI. The telescope 
was defocussed slightly to keep the counts for all the stars 



Instrument 


Epoch a 


T c (BJD TDB ) b 


O - C (s) c 


MMIRS 


464 


2455699.832866 ± 6.9E-5 


-7.6±5.9 


FORS 


493 


2455745.664729 ± 6.0E-5 


+3.1 ±5.2 


HAWKI 


517 


2455783.59440 ±l.lE-4 


-0.4±9.5 



a Integer number of orbital periods since the transit on 2009 May 
15. 

BJDtdb is the the Barycentric Julian Date in the Barycentric 
Dynamical Time standard (Eastman ct al.ll2"010l ). 
c Residuals form the ephemeris given in §3. 



in the field-of-view within the linearity range of the de- 
tector. To improve the duty cycle, six separate exposures 
were averaged on the chip to produce a single recorded 
image. The total time to integrate the exposures, read 
out the data, and reset the detector was 28 s. A total 
of 337 images were obtained, 110 of which were during 
transit. The maximum counts for GJ 1214, which was 
the brightest star in the field-of-view, ranged from 9,000 
to 25,000 e~. This is well below the the 1% non-linearity 
level for the HAWKI detectors (approximately 60,000 e~) 
and no non-linearity corrections were applied to the data. 
The seeing reported by the observatory atmospheric sta- 
tion varied from 0.9 to 3.2" during the observations, but 
we recall that the guide probe of the telescope never reg- 
istered more than 1.6" (this information is not included 
in the image headers). 

After bias subtraction and flat fielding the HAWKI 
images, we measured aperture photometry for GJ 1214 
and 11 other stars. The photometric aperture and sky 
annulus sizes were optimized to give the lowest dispersion 
in the light curve for GJ 1214. We settled on a stellar 
aperture radius of 26 pixels, and a sky annulus inner 
radius of 31 pixels and outer radius of 41 pixels for the 
final analysis. 

We corrected the time-series photometry of GJ 1214 
by dividing out the sum of a number of the reference 
stars. We experimented with different combinations of 
the reference stars for the correction, and we ultimately 
used a composite of four of them for the final analysis. A 
plot of the raw time-series for GJ 1214 and the composite 
reference is shown in Figure [5] 

After the reference star correction, we noticed a slow 
and smooth variation in the corrected GJ 1214 relative 
flux values with time superimposed on the expected tran- 
sit light curve morphology. This trend is potentially re- 
lated to the large color difference between GJ 1214 and 
the reference stars, although the bandpass of the utilized 
filter is relatively narrow. Similar trends are regularly 
seen in ground-based n ea r-infra re d transit photom etry 
(e.g.. ICroll et~alll2010allR I2011allb1 : ISada et alJl2010h . so 
this probably arises from a property of Earth's atmo- 
sphere. We found that the trend in our data was best 
modeled by a second order polynomial with time, and 
we fit for such a function simultaneously with the transit 
modeling (see §3, two free parameters in addition to the 
normalization). Similar results, but with higher disper- 
sion in the light curve residuals and larger determined 
parameter uncertainties, are obtained when decorrelat- 
ing the data with a linear trend as a function of airmass. 
The data from three of the images were found to deviate 
by more than five times the standard deviation of the 
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Fig. 7. — Upper panels Normalized broadband (left panel) and TC-band spectroscopic (right panel) light curves from the MMIRS data 
(circles). The best-fit models are shown as the solid lines. Bottom panels Residuals from the fits (circles). 



residuals from a best-fit model. These discrepant points 
were removed for the final analysis. 

The HAWKI light curve for GJ 1214 after all correc- 
tions, decorrelations, and trimming of discrepant points 
exhibits residuals from a best-fit model with a rms of 
1259 ppm. The rms of the residuals is 732 ppm when the 
data are binned to 60 s. The dispersion of the per-image 
residuals in the un-binned data is a factor of 1.3 larger 
than expected from the estimated photon noise. 

3. ANALYSIS 

We fit the broadband and spectroscopic light curves 
for GJ 1214 with a transit model multiplied by a nor- 
malization factor and systematic decorrelation functions 
to investigate the transmission spectrum of its transiting 
planet. We used the exact anal ytic formulas including 
quadratic limb darkening given by Mand el fc Agoll (2002) 
for the transit model. The transit model was parameter- 
ized by the planet-to-star radius ratio (Rp/R*), system 
scale (a/R-t), orbital inclination of the planet (i), mid- 
transit time of the planet (T c ), and the stellar quadratic 
limb darkening coefficients (71 and 72). The orbital pe- 
riod of the planet was fixed to the revised value given 
below, and the eccentricity was fixed to zero. The nor- 
malization factor for each data set was always a free pa- 
rameter. The decorrelation functions described in §2 re- 
quire two additional free parameters for each wavelength 
bin. 

We estimated limb darkening coefficients appropriate 



for transits of GJ 1214b using spherically symmetric stel- 
lar m odel atmospheres calcu lated with the PHOENIX 
code (jHauschildt et al.lll999| ). We adopted the stellar 
parameters for the baseline model T e ff = 3026 K, [M/H] 
= 0.0, logg = 5.0, and R+ = 0.21 i? Q . We also computed 
models at the baseline effective temperature ± 130 K. For 
the MMIRS and FORS spectroscopic data, the model 
intensities at each wavelength were scaled to the actual 
counts in the utilized spectra for GJ 1214 to account for 
the unique bandpass and flux weighting of the resulting 
photometry. For the HAWKI photometric data, we es- 
timated the effective observed stellar spectrum by mul- 
tiplying the baseline theoretical stellar spectrum by a 
theoretical model for the transmission of Earth's atmo- 
sphere and by the transmission curve of the utilized fil- 
ter. The theoretical model for the transmission of Earth's 
atmosphere was calculated for the line of sight through 
atmosphere above the observatory on the night of the ob- 
servat ions using the method described bv lSeifahrt et all 
(2010). After integrating over the effective bandpasses, 
we fit a quadratic function to the intensities to determine 
the limb darkening coefficients. We limited our fits to the 
range /i = cos 9 < 0.1, where 6 is the angle between the 
emergent intensity and the line of sight. 

We allowed the limb darkening coefficients to be free 
parameters in the light curve modeling, and we used 
the estimated theoretical values as priors. That is, the 
differences between the coefficients used for the light- 
curve modeling and the coefficients calculated from the 
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TABLE 3 

Photometric Bandpasses & Transit Depths 
for the FORS Blue Data 



Band 



Wavelength (nm) 



Fig. 8. — Left panel Normalized light curves from the FORS blue data (circles) and best-fit models (lines), 
the fits (circles). 

PHOENIX models were included in the tabulation of the 
goodness-of-fit metric (see below). The adopted uncer- 
tainties for the priors were set by the spread in the values 
for the coefficients at the three different temperatures 
considere d for GJ 1214 . We also tested the effects of 
using the IClaretl C2000) non-linear limb darkening law. 
While this functional form yielded better fits at fi < 0.1 
for the MMIRS data, it did not yield better light curve 
fits and the final determined transmision spectrum val- 
ues for the planet are not significantly affected by using it 
instead of the quadratic law. The best-fit limb darkening 
coefficients were in all cases within 2 a of the theoretical 
values estimated from the T e // = 3026 K model atmo- 
sphere. 

We first analyzed the broadband light curves (MMIRS 
J, H, and K- FORS blue and red; and HAWKI) to de- 
termine values for all the transit parameters. We ana- 
lyzed the MMIRS, FORS, and HAWKI data separately. 
We assumed the system scale, inclination, and transit 
time are the same for the J-, H-, and if-band MMIRS 
data because they were obtained simultaneously for the 
same transit. Each of the MMIRS bands were allowed to 
have a unique tran sit depth. We use d a non-linear least- 
squares algorithm (Markwardt 2009) to identify the best- 
fit parameters, and a residual permutation bootstrap al- 
gorithm to asses the uncertainties on the parameters. 
The standard \ 2 metric was used throughout to asses 
the quality of the model fits. 

The determined transit times from analysis of the 



Rp/ R+ 



broadband 


610 


850 


0.1178 ±0.0007 


channel 1 


610 


630 


0.1173 ±0.0018 


channel 2 


630 


650 


0.1195 ±0.0012 


channel 3 


650 


670 


0.1167 ±0.0011 


channel 4 


670 


690 


0.1194 ±0.0011 


channel 5 


690 


710 


0.1169 ±0.0009 


channel 6 


710 


730 


0.1164 ±0.0009 


channel 7 


730 


750 


0.1182 ±0.0007 


channel 8 


750 


770 


0.1187 ±0.0008 


channel 9 


770 


790 


0.1172 ±0.0008 


channel 10 


790 


810 


0.1172 ±0.0007 


channel 11 


810 


830 


0.1183 ±0.0006 


channel 12 


830 


850 


0.1168 ±0.0007 



a From an analysis with the transit parameters a/R* 
and i fixed to 14.97 and 88.94°, respectively. 



broadband data with all parameters free are given in 
Table [U Combinin g these times with t he previously 
report ed ti mes from iDesert et al.l (|2011l) . iCarter et al.l 
(pOll . and IBerta et al.l (|2011D . we can calculate a re- 
vised period P = 1.58040481 ± 1.2E-7d, and reference 
transit time T c = 2454966.525123± 0.000032 BJD T db- 
The new transit times in this paper deviate by less than 
1.3t7 from the predictions of this ephemeris, and none of 
other transit times exhibit deviations greater than 2.4 a. 
The best-fit values of the system scale and the planet's 
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TABLE 4 

Photometric Bandpasses & Transit Depths 
for the FORS Red Data 



Band 



Wavelength (nm) 



Rp/ Ft* 



broad band 


7so 


i ooo 

1UUU 


11 fiS 
U. 1 lOO 


ztz 


O OOOfl 
U.UUUD 


channel 1 


7£0 


7Q0 


nil fi7 


4- 
zrz 


n oooq 
u.uuuy 


channel 2 


7Q0 
/ yu 


Ron 

OUU 


U. 1 1UU 


4- 

ZL 


n ono7 

U.UUU ( 


channel 3 


son 

OUU 


si o 

OlU 


U. 1 lOO 


4- 

ztz 


O 0007 
U.UUU i 


channel 4 


sin 

OlU 


S90 


11 7(\ 
U. 1 1 ( D 


4_ 

ZtZ 


ooos 
u.uuuo 


channel 5 


&90 
ozu 
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1 1 7fi 
U. 1 1 ID 




0007 
U.UUU 1 


channel 6 


OOU 


SAO 

o^tu 


11 fi9 
U. 11DZ 


ztz 


O 0007 
U.UUU i 


channel 7 


S40 

O^tU 


s^o 

OOU 


11 79 


_l_ 

ZtZ 


O 0007 
U.UUU i 


channel 8 


850 


860 


0.1151 


± 


0.0008 


channel 9 


860 


870 


0.1168 


± 


0.0007 


channel 10 


870 


880 


0.1171 


± 


0.0007 


channel 11 


880 


890 


0.1171 


± 


0.0007 


channel 12 


890 


900 


0.1159 


± 


0.0007 


channel 13 


900 


910 


0.1167 


± 


0.0006 


channel 14 


910 


920 


0.1175 


± 


0.0007 


channel 15 


920 


930 


0.1178 


± 


0.0006 


channel 16 


930 


940 


0.1165 


± 


0.0009 


channel 17 


940 


950 


0.1168 


± 


0.0008 


channel 18 


950 


960 


0.1176 


± 


0.0009 


channel 19 


960 


970 


0.1172 


± 


0.0009 


channel 20 


970 


980 


0.1161 


± 


0.0010 


channel 21 


980 


990 


0.1165 


± 


0.0011 


channel 22 


990 


1000 


0.1168 


± 


0.0011 



a From an analysis with the transit parameters a/R* 
and i fixed to 14.97 and 88.94°, respectively. 

orbital inclination for the broadband data are consistent 
with previously determ ined values (e.g., iCarter et al.1 
120111 iBerta et al.ll2011|). We fixed the se parameters to 
the values used by iBean et all (|201Ct a/R* = 14.9749, 
i = 88.94°) to measure the transit depth as a function 
of wavelength in the subsequent analyse s. T hese values 
were al so adopted bv lDesert et al.l (|2011l ) and lCroll et al.l 
(|2011al ) in their studies, so we can directly compare our 
results with theirs. We also fixed all the transit times 
to the predicted values (not best-fit) from the ephemeris 
given above. This is justified because there is no evidence 
of transit timing variations for this system. 

We performed two types of analyses to measure the 
planet's transmission spectrum. We used non- linear 
least-squares fitting with a residual permutation boot- 
strap to analyze the broadband light curves. For the 
spectroscopic data sets, we used a Markov Chain Monte 
Carlo (MCMC) algorithm to find the most likely param- 
eter values and their uncertainties. The adopted photo- 
metric errors for the MCMC analysis were the photon- 
limited errors adjusted upwards to yield a reduced \ 2 = 1 
for the best fit. The free parameters for each wavelength 
bin in the two analyses were R p /R+, quadratic limb- 
darkening coefficients (with priors), flux normalization, 
and decorrelation parameters. There were a total of six 
free parameters for each wavelength bin, and two "obser- 
vations" per bin in addition to the light curves (the limb 
darkening priors). 

The motivation for the different analyses between the 
broadband and spectroscopic data sets is that they are 
used in different ways to investigate the planet's trans- 
mission spectrum. The broadband data are useful when 
compared to theoretical models in aggregate because 
they span a wide range of wavelengths. For this kind 
of study, the correlated, or red, noise in each data set 
must be considered because the observations were ob- 
tained for different transits with different instruments 
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Fig. 9. — Upper panel The normalized HAWKI light curve (cir- 
cles) and best-fit model (line). Lower panel Residuals from the 
best-fit model (circles). 

over the course of years. On the other hand, the spectro- 
scopic data sets each tightly constrain the planet's trans- 
mission spectrum on their own. These measurements 
were obtained simultaneously with the same instrument 
for the same transit and over a limited wavelength range. 
In this case, the major sources of correlated noise, vari- 
ations in Earth's atmospheric transparency and stellar 
spot crossing, arc quite similar from channel to channel 
within a data set. Therefore, a full correlated noise anal- 
ysis is not appropriate when the resulting planet trans- 
mission spectra for each spectroscopic data set are only 
examined in isolation, as we do in §4. 

The determined R p /R± values for the MMIRS, FORS 
blue, and FORS red data are given in Tables [TJ [3J and 
[31 respectively. The determined R p /R+ for the HAWKI 
data is 0.1179 ± 0.0012. The normalized light curves with 
systematics removed along with the best-fit transit mod- 
els and residuals from the models are shown in Figures [7] 
(MMIRS data), ® (FORS blue data), and M (HAWKI 
data). The re-analyzed FORS red data do n ot look sub- 
stanti ally different than those data shown by IBean et all 
(|201d spectroscopic data) and IBerta et al.l ( 20111 . broad- 
band data). 

There is some indication of the planet crossing spots on 
the surface of GJ 1214 in the FORS blue data. We per- 
formed tests to investigate the influence of masking the 
possible spot crossing events from the light curves. The 
determined depths of the light curves changed depend- 
ing on which points were masked, but the variation was 
always less than the quoted uncertainties on the R p / R+ 
values, and the relative values of the FORS blue trans- 
mission spectrum were unchanged. The results presented 
here are from an analysis with none of the light curve 
points masked. 

4. DISCUSSION 

We compare our newly measured transit depths to 
previous results and the predictions of theoretical mod- 
els to investigate the nature of GJ 1214b's atmosphere 
in this section. We u tilize the models presented by 
iMiller-Ricci fc Fortnevl (|2010T > throughout. The goals of 
our investigation are to determine the composition of the 
planet's atmosphere, and to provide an independent as- 
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sessment o f the transit depth in the -ftT-band to test the 
claim from ICroll et all (|2011aD that the transit is signif- 
icantly deeper at these wavelengths. 

We look at the broadband, and the optical and K- 
band spectroscopic data sets separately. The spectro- 
scopic data enable a search for spectral features within 
the different regions they cover. The strength of these 
data sets is that they yield extremely precise measure- 
ments on an internal scale, and each one tightly con- 
strains the planet's transmission spectrum. However, 
there are offsets between the different spectroscopic data 
sets with sizes that are on the order of the uncertain- 
ties of the respective broadband points. The FORS blue 
and red data can be combined because the offset between 
the data sets can be determined from the values at the 
wavelengths where they overlap. The FORS and MMIRS 
spectroscopy can not be combined because the offset be- 
tween the data sets is uncertain. 

The broadband data cover a wide span of wavelengths, 
and thus probe different sources of opacity and dif- 
ferent atmospheric pressures. These data were taken 
at different epochs and with different instruments, but 
can be combined because we considered the correlated 
noise in each of the data sets, and because the activity- 
related photometric variability of GJ 1214 itself should 
only cause variations in the transit depth th at are much 
small er than the measurement uncertainties (jBerta et al.1 
2011). The excellent agreement between the FORS blue 
and red observations, which were separated by a year in 
time, at the wavelengths the data sets have in common 
(see §4.3) is additional evidence that stellar activity is 
not a significant problem for this investigation. 

4.1. Broadband photometry 

The broadband m easurements o f GJ 12 14b's transmis- 
sion spectrum fro m IDesert et al.1 (|201lL Spitzer data) , 
ICroll et ahl (|2011al ground-based near-infrared), and this 
paper are shown compared to theoretical models for the 
planet's atmosphere in Figure [TO] The broadband trans- 
mission spectrum from the combination of the FORS 



blue and red, MMIRS, HAWKI, and Spitzer data does 
not exhibit any significant features. The maximum devi- 
ation from the weighted mean of the points is 1.7 a. 

Up until now, the only candidate detection of spectral 
features in GJ 1 214b's transmission spectrum has been in 
the if-band by ICroll et aD (|2011aD . Our HAWKI data 
point in a narrow window of the K-band do es not have 
the pre cision to provide a strong test of the ICroll et "all 
(2011a) measurement on its own. However, the broad 
-ftT-band data point created from the combination of the 
MMIRS spectroscopic data has excellent precision and 
can provide this test. The MMIRS if-band p oint is 4.1 a 
A R„/ R* = 0.0041 ± 0.0010) lower than the ICroll et "all 
2011af) if-band point, and is consistent with all of our 
group's measurements at other wavelengths. Therefore, 
we do not c onfirm the detection of spectral features by 
ICroll et all (|2011aD . 

The theoretical models for the planet's atmosphere 
we compare to the observations were scaled to give the 
best-fit to the FORS blue and red, MMIRS, HAWKI, 
and Spitzer points. There are eight data points and 
one free parameter, and thus seven degrees of freedom 
for the model fits. A model for the planet's atmo- 
sphere with a 100% water composition, which is essen- 
tially fiat at the level of the precision of the measure- 
ments, yields x 2 = 10.1. On the other hand, a model 
for the planet's atmosphere with a solar composition 
(i.e., hydrogen-dominated) and assuming chemical equi- 
librium gives x 2 = 83.9. This suggests that the simple 
solar composition model is ruled out at 7.9 a confidence. 
This strengthens the conclusions we r eached previously 
using subsets of t he current data set (|Bean et al.l 12010 : 
IDesert et al1l20Tfl) . 

Methane is an important opacity source for solar com- 
position models of GJ 1214b's cool atmosphere, partic- 
ularly in the if-band and at 3.6 /xm. However, the 
predicted abundance of this molecule is subject to sig- 
nificant uncertainty because it could be affected by 
non-equilibrium processes (e.g., photochemistry, ther- 
mal chemistry, and mixing). Millcr-R icci Kempton et al.l 
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(|2011[ ) studied this possibility and concluded that 
methane depletion can at most only reach atmospheric 
pressure levels of 10 -4 bar, which is far above the al- 
titude probed by transmission spectroscopy (approxi- 
mately 0.01 to lbar). Similar conclusions have been 
reached about the difficulty of photochemistry to sig- 
nificantly alter the methane abundance at observable 
depths in the atm osphere of the warmer planet GJ 436b 
(jLine et al.ll2011| ) . Comparing a solar-composition model 
with methane artificially removed to our group's broad- 
band data for GJ 1214b gives x 2 — 43.8, which is dis- 
crepant at 5.2 a. 

A grey broadband optical opacity source is needed for 
hydrogen-dominated models of the planet's atmosphere 
to be consistent with the data. Such a source could be 
high-altitude clouds or haze. The data constrain the al- 
titude of clouds or haze in GJ 1214b's atmosphere to be 
at least above a level corresponding to pressures of ap- 
proximately 0.1 bar. If methane is depleted down to a 
similar level, then the clouds or haze would only have to 
be optically thick for wavelengths less than 1 /im. 

4.2. K-band spectroscopy 

In addition to providing broadband measurements, the 
MMIRS data also yield a low-resolution spectrum of 
GJ 1214b in the if-band. These data are shown com- 
pared to theoretical models in Figure QTJ We examine 
these data in isolation, and adjust the models to give the 
best fit without consideration for the measurements at 
other wavelengths. There are three degrees of freedom 
for the model comparisons. 

The MMIRS if-band spectrum is featureless and fiat 
like the broadband data. The maximum deviation from 
the weighted mean of the points is 0.5 a. So not only 
do we not confirm in our b roadband data the d eeper K- 
band transit suggested by iCroll et al.l (|2011aD . we also 
do not detect any spectral features in the planet's at- 
mosphere at_Jh^se__wavdeii£ths using differential spec- 
troscopy. iCrossfield et al.l ()2011l ) previously presented a 
non-detection of spectral features in this region based on 
high-resolution differential spectroscopy, and we confirm 
their results at lower resolution. 

The 100% water atmosphere model is consistent with 



the data (x 2 — 0.9), while the solar composition model in 
equilibrium is ruled out at 4.3 a confidence using these 
data alone (x 2 — 24.5). The main opacity source in this 
window is methane, and the solar composition model 
with methane artificially removed is consistent with the 
data (x 2 — 2-0). The conclusion we draw from these data 
is that the planet's atmosphere must either be metal- 
enhanced or mcthanc-dcplctcd. 

4.3. Optical spectroscopy 

The optical spectroscopic measurements from the new 
FORS blue data and the re-analyzed FORS red data are 
shown with theoretical models in Figure [T2"1 We examine 
these data in isolation, and adjust the models to give the 
best fit without consideration for the measurements at 
other wavelengths. We find that the FORS blue values 
are on average 0.0007 larger (indicating a deeper transit) 
than the FORS red data for the wavelength range the 
two data sets have in common. This is less than the 
lcr uncertainty (0.0011) between the broadband points 
for each data set. We subtracted this determined offset 
from the FORS blue data for display in Figure [12] and 
the calculation of the fit quality for the different models. 
The FORS optical data include 34 measurements, and 
there are 33 degrees of freedom for the examination of 
the model fit quality. 

The maximum deviation from the weighted mean of 
the FORS spectroscopic points is 2.3 a, which suggests 
that no spectral features are detected in these data. 
The solar composition model in equilibrium is ruled 
out at 6.6 a confidence (x 2 = 115.7) using these data 
alone. The solar composition model with methane re- 
moved is discrepant from the FORS data at the 5.2 a 
level (x 2 = 91.2). The 100% water composition model is 
consistent with the data (x 2 = 32.4). Water mass frac- 
tions of 70% are needed to bring the models within 3 a 
of the data, suggesting a highly metal-rich atmosphere. 
This conclusion is the same as we found b efore using 
only the FORS red data at lower resolution (tBean et al.l 
I2010T ) . The higher resolution and precision for the FORS 
red data enabled by our improved data reduction, and 
the addition of the blue data tightens the constraints on 
the planet's optical transmission spectrum. 

5. CONCLUSIONS 

We have obtained new ground-based measurements of 
the transmission spectrum of the 6.5 M® planet GJ 1214b 
in the optical between 0.61 and 0.85 fim, and in the J, H , 
and K near-infrared windows. We have also re- analyzed 
our previously reported red optical spectroscopy for the 
planet between 0.78 and 1.00 /xm. We were able to 
push to higher resolution transmission spectrum mea- 
surements with these data because of improvements in 
the data reduction algorithm. 

We combined the new data with previously reported 
measurements for GJ 1214b's t ransmission spectru m that 
were obtained with Spitzer (|Desert et alj 120 III) . The 
combined data set spans the visible to the infrared 
(0.6 to 4.5 /jm), which makes it one of the most com- 
plete exoplanet transmission spectra obtained to date. 
We compared the combine d data set t o the p reviously 
reported measurements of Cro ll et all ()2011al ) and to 
theoretical models for the p lanet's atmosphere from 
iMiller-Ricci fc Fortnevl (|2010j ). Our main conclusion is 
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that there is no evidence of features in the planet's trans- 
mission spectrum, and we do not confirm the detec- 
tion of a significan tly deeper transit in the if-band by 
ICroll et aLll20TTal) . 

Our current knowledge of GJ 1214b's mass and ra- 
dius indicates that the p l anet must have an atmo sphere 
(|Rogers fc Seagerj|2010bT ; INettelmann et al.ll201lD . Our 
interpretation of the featureless transmission spectrum 
for the planet is that its atmosphere must either have at 
least 70% H2O by mass or optically thick high-altitude 
clouds or haze. A simple cloud- and haze-free model for 
the planet's atmosphere with solar composition gas in 
chemical equilibrium is ruled out at high confidence. 

Alternatively, our knowledge of the planet's radius 
could be wrong. If the planet were actually significantly 
smaller (approximately 2R^, or 6.5 a from the current 
best estimate), then it would not necessarily need to have 
a substantial atmosphere. The estimate of GJ 1214b's 
radius depends on the assumed distance to the system, 
which is used to estimate the mass of GJ 1214 itself, and 
this is the link in the chain with the biggest question 
mark above it because the most recent estimate is based 
on photographic plate measurements (|van Altena et all 



1995). If the parallax were approximately 131 mas (10 a 
greater than the current estimate), the radius of the 
planet determined from the light curves would be con- 
sistent with the planet having no atmosphere. A new 
measurement of the system's trigonometric parallax us- 
ing modern technology would further our knowledge of 
this important planet. 
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